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Abstract 



We summarize studies of growing lengths in different aging systems. The article is 
structured as follows. We recall the definition of a number of observables, typically 
correlations and susceptibilities, that give access to dynamic and static correlation 
lengths. We use a growing length perspective to review three out of equilibrium cases: 
domain growth phenomena; the evolution of Edwards- Wilkinson and Kardar-Parisi- 
Zhang manifolds and other directed elastic manifolds in random media; spin and 
structural glasses in relaxation and under an external drive. Finally, we briefly re- 
port on a mechanism for dynamic fluctuations in aging systems that is based on a 
time-reparametrization invariance scenario and may be at the origin of the dynamic 
growing length in glassy materials. 



1 Introduction 



For a long time, and somewhat paradoxically, the majority of theoretical physicists 
interested in glasses were reluctant to study these systems in the trully glassy regime, 
say below the glass temperature T g or above the glass density p g . One of the reasons 
to resist entering the glassy regime was the lack of insight on which questions to ask 
and, more concretely, which quantities to measure in out of equilibrium conditions. 

The situation changed dramatically around 20 years ago when it was accepted that 
glassy systems do not freeze out below T g but they just continue to relax at a slower 
and slower rate as time goes by - the aging phenomenon. Although this fact was 
well-known experimentally and to a certain extent also numerically it was not fully 
assimilated theoretically. The solution to a number of simple models [namely, droplet 
models for disordered systems (JJ), random walks among traps in phase space ([2]), 
and mean-field spin systems with disordered interactions ([3] SJ] demonstrated the 
possibility of capturing aging phenomena analytically. These solutions also provided 
a guideline as to which are the simplest preparation protocol and the most relevant 
obscrvables to focus on dynamically. More importantly, these studies also opened the 
way to confront glassy dynamics to the behaviour of other macroscopic non-equilibrium 
systems considered to be simpler such as coarsening phenomena ([5]) or the motion of 
manifolds in different kinds of embedding volumes ©• Experimental and numerical 
data on theoretically motivated measurements in various physical systems have rapidly 
accumulated and some of these results are summarized in other chapters in this book. 

As for the measuring protocols, in short, one first chooses the way in which the 
system is taken into the glassy regime and defines the origin of the time axis accord- 
ingly. In the simplest setting, the 'initial time' t = 0, is the instant when the system 
is suddenly prepared in an out-of-equilibrium condition. For instance, this is achieved 
by performing a rapid temperature quench from a high temperature above T a down 
to a target value T <T g but other routes to the glassy regime, changing other control 
parameters, are also feasible. As will be discussed below, the physical properties of the 
system become functions of the waiting time t w , the time spent after the preparation 
of the initial state, stationarity is lost and the relaxation time increases with t w . The 
basic idea to analyze such an aging relaxation is to introduce a laboratory time scale 
- the waiting-time - that can be controlled at will when the intrinsic time scale - the 
equilibrium relaxation time - goes beyond accessible times. 

Following this procedure one can study, for instance, the energy relaxation as a 
function of t w . However, such a macroscopic 'one-time quantity' does not capture all 
the richness of the relaxation. Much more detailed information is contained in two 
(or more) time-dependent global quantities. In Fig. ll.ll we show the typical behaviour 
of two macroscopic two-time quantities in many glassy systems: the spontaneous de- 
correlation of a chosen global observable measured at two subsequent times t and t w 
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Fig. 1.1 Typical behaviour of macroscopic two-time quantities during aging, (a) Auto-corre- 
lation function C(t,t w ); (b) linear-susceptibility \(t,t w ). The waiting time t w increases from 
the left to the right curves, typically in a logarithmic scale such that, say, t W k+i = 10£ m fc. 
The a relaxation time, t a , or the time-difference needed to decorrelate significantly, is an 
increasing function of t w and is indicated with an arrow in the left panel. The plateau in \ 
occurs at \ ea = (1 — q? a )/T where q ea is the height of the plateau in C. 

(a), and the linear response of the same observable measured at t to a perturbation 
that couples linearly to it between t w and t (b). 

As shown in panel (a) , the auto-correlation function exhibits a two-step decay. At 
short time scales, say t — t w <C t w , there is a relaxation towards a plateau that we call 
the Edwards- Anderson (EA) order-parameter q ea (see its precise definition in Eq. (O). 
The relaxation in this stage is essentially independent of the waiting time t w and time- 
translational invariance (TTI) holds C(t,t w ) = C s t(t — t w ) to a certain accuracy. The 
second relaxation at longer time scales, say t — t w ^> t w , exhibits i^-dependence and 
is not stationary. The latter reflects the strong out-of equilibrium nature of aging. The 
separation of time-scales is confirmed in panel (b) where we sketch the behaviour of 
the susceptibility x{t^w)- The plateau occurs here at Xea = (1 — Qea)/T with q ea the 
value of C at its plateau and T the temperature of the environment that is henceforth 
measured in units of fee, the Boltzmann constant. 

The relaxation of the correlation in Fig. 11.11 (a) reminds us, in a sense, of the 
two step relaxation in the supercooled liquid phase - the so-called a and f3 relax- 
ations ([7]). Even more so, it is basically indistinguishable from the one found in phase 
ordering processes after a temperature quench (see Sect. [3] for a detailed descrip- 
tion of coarsening) ([5]). In such processes domains of the different equilibrium phases 
progressively grow in competition. At each instant, the equilibrium order parameter is 
essentially uniform within each domain. When only short time-differences are explored 
the domain-walls are basically static and the correlations decay just as in equilibrium. 
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Instead, when longer time- differences are reached the walls move appreciably and the 
correlations decay from the plateau in a manner that depends explicitly on the waiting- 
time. A key quantity in the description of coarsening is the typical domain linear size, 
L(t), which plays the role of a dictionary between time t and length L. All properties of 
coarsening systems are invariant under rescaling of all lengths by L(t) - in a statistical 
sense. 

The question naturally arises as to whether the dynamics in glassy systems is also 
of coarsening type albeit the growing order had not been identified yet. From the 
analysis of the correlations in many glasses, and one example is shown in Fig. 11.21 
where scaled data from an aging colloidal suspension are displayed (O, one could 
easily conclude this to be the case. The figure shows the two-time correlation data 
taken using different t w s as a function of the ratio t/t w . For sufficiently long t w the 
scaling is quite good suggesting that there might be a growing length Lit) oc t a , with 
a an undetermined power. However, the scaling of the global two-time correlation 
as f (L(t) / L(t w )) is just a quite generic property of two-time correlations (j4]), and 
the space-time correlations do not signal the growth of structure contrary to what 
happens in simple phase ordering kinetics in which it scales as a function of r/L(t). 
Therefore, to reach a conclusion on the existence (or not) of a growing length one 
should, ideally, analyze the dynamics at mesoscopic time and length scales or, at least, 
define and explore more complex quantities giving access to the mesoscopic details of 
the process. Eventually, these studies might also enlighten us about the reason for the 
growth. 




t/t w 



Fig. 1.2 (a) Scaling of two-time correlation data in an aging colloidal suspensions. Waiting 
times are given in the key in units of 10 2 s. The relatively good scaling with t/t w suggests 
L(t) ~ t a . The data are taken from and analyzed along the lines explained in ([SJ. (b) A 
sketch of the pixelization used to describe the particle configurations with binary variables, 
a a = 0, 1; and the coarse-graining procedure (see Sec. 12. a is the lattice spacing, D the 
particle diameter, £ the coarse-graining length and L the system size. 
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Taking the case of usual phase ordering processes as a reference, we realize that 
in order to settle this issue we need to follow a number of basic steps summarized 
by the following questions: 1) how can wc define an order parameter for a glass? 2) 
how can we define domains - within which the order parameter is essentially constant 
- or dynamically correlated volumes? 3) can one identify a growing linear size and 
does it reflect the nature of the dynamics - e.g. smooth (such as, for instance, the 
curvature driven kinetics observed in some coarsening systems) vs. activated? 4) can 
one understand the mechanism leading to such a growth? In the body of this chapter 
we develop these questions. 

Having explained that the aim of this chapter is to discuss searches of growing 
lengths in aging glasses, we announce that we shall also deal with a priori simpler cases 
such as a particle's diffusion in random media or the motion of a finite dimensional 
manifold, both in their aging regime. Although the glassy nature of these problems 
might not be complete, it turns out that some aspects of their dynamics are similar 
to the ones encountered in conventional glasses. Moreover, their mere definition is 
clearer and their analytic treatment as well as the identification of a growing length 
are simpler. Their discussion should then be instructive. 

The organization of the chapter is the following. In Sect. [5] we give a number of 
definitions and we discuss how to deal with quenched disorder, if present. In Sect. [3] 
we review sub-critical and critical coarsening phenomena in clean systems. In this 
Section we also introduce, briefly, a relatively simpler problem with an easy to identify 
growing length: the motion of a free elastic manifold. Section 0] is devoted to the 
phcnomenological decription of activated processes - that are specially relevant in 
problems with quenched disorder - starting with a short account of aging in the Sinai 
model of difussion and developing afterwards the droplet model predictions for the 
dynamics of clastic manifolds in random media and spin-glasses. Section [5] recaps 
mean-field predictions for the growing length issued from the study of fully-conncctcd 
spin models and the related mode-coupling theories as well as growth phenomena in 
kinetically constrained models for glasses. We also discuss growing lengths in other 
systems with aging dynamics: the interplay between flow and internal relaxation and 
the effect of the former on the aging properties of coarsening and glassy systems, 
granular matter and quantum glasses. Finally, in Sect. [6] wc briefly recall a possible 
mechanism for dynamic fluctuations and the associated growing correlation length in 
aging systems: the asymptotic development of time-reparametrization invariance, that 
has been extensively reviewed elsewhere ([9]). 

We do not attempt to present a comprehensive review of analytic, numeric and 
experimental studies of all kinds of aging systems. We rather focus on the search for 
growing length-scales and their possible use to describe aging materials. With this 
in mind, the reader should not expect to find a complete list of references on aging 
studies. 

2 Definitions 

In this Section we summarize a number of definitions relevant to the study of cooper- 
ative motion in glassy systems. These definitions do not assume equilibrium and can 
be used to study out of equilibrium and, in particular, aging samples. They depend, 
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in general, on various times independently and no fluctuation-dissipation theorem is 
assumed. Although we are not going to dwell into the question of the pertinence of 
effective temperatures in the description of glasses, we recall its definition since it will 
appear in the rest of this Chapter. 

2.1 Euler and Lagrange descriptions 

Two ways of studying the heterogeneous dynamics in a many particle (or many higher- 
dimensional object) system are the following. 

Lagrangian description. One can follow the evolution of each individual particle, la- 
beled by an index i, and detect which are the fast and slow moving ones during a 
previously chosen time window, say t — t w , around some time, say t w , after prepara- 
tion. This is the route followed in early studies of dynamic heterogeneities in super- 
cooled liquids, see (|10p and (fTTj) for a few of many molecular dynamics and confocal 
microscopy papers that use this type of measurement. One of the main outcomes of 
these studies is the observation of clustering of fast particles, in the form of strings the 
length of which increases - and may diverge - close to the glassy arrest. In the trully 
glassy phase fewer studies exist. A precursor is (fT2j). More recently, Vollmayr-Lee et al 
studied the geometry and statistical properties of clusters of mobile and inmobile par- 
ticles with molecular dynamics (13). Confocal microscopy has been used with similar 
aims CHI)- 

Eulcrian description. In a first step, one can define Ising spin variables as the ones 
that describe uniaxial magnetic systems. The choice is done for notation convenience 
but also because, as proposed in ©, a simple mapping between particle positions 
and Ising spin variables, o~ a , defined on the N = (L/a) 3 vertices of a cubic lattice 
captures the dynamics of 'atomic' glassy systems as well (see Fig. 11.21 (b)). (L is the 
linear length of the box and a the length of the lattice spacing.) In a few words, one 
partitions space with a lattice with very fine mesh - smaller than the particle radii - 
and the mapping assigns a spin one to each cell occupied by a piece of particle and 
zero otherwise. Once this is done, all observables are written in terms of the bimodal 
variables, as in spin models, and the identities of the particles are ignored. 

In a second step, one can coarse-grain the spin variables over boxes of a chosen 
linear size, £, with ideally L ;§> £ ;§> a. The locality is then given by the position of the 
box which is labeled by an index i. The coarse-graining amounts to a partial averaging, 
the effective spin on each box, Si, is now a continuous variable, < s, < 1. 

In a system with quenched randomness, realized as fixed obstacles, local fields 
or else, the coarse-graining procedure over a sufficiently large volume averages over 
the peculiar local features of the fixed disorder. In this way, the 'finger-print' (|15p of 
disorder, that is to say local fluctuations determined by the particular local disorder - 
such as Griffiths singularities in random ferromagnets - should be washed away. One 
can then expect to arrive at a description of the noise-induced fluctuations present in 
problems with or without quenched disorder. 

This method appears to be more adequate for analytic treatment through a field 
theory. 
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2.2 Space-time correlation 

In usual coarsening systems, see Sect. El the averaged space-time correlation function 
NC(r,t)= Y, (Ssi(t)5 Sj (t)} , (1) 

ij/\ri-rj\=r 

with 8si(t) = Si(t) — (si(t)) and (si(t)) = (si) eq in all the cases we shall deal with, allows 
for the identification of a growing length from, for example, L a (t) = J d d r r a+1 C(r, t) / 
J d d r r a C(r, t). (a is a parameter chosen to weight preferentially short or long distances; 
the time-dependence of L a (t) should not depend on a.) Here and in the following 
(. . .) stands for an average over different realizations of thermal histories at heat- 
bath temperature T and/or initial conditions. In presence of quenched disorder one 
adds an average over it and denotes it [...]. The stochastic time-dependent function 
J2ij/\ri-r \=r s i(^) s j(^) a fter a quench from a random initial condition does not 
fluctuate in the thermodynamic limit. Therefore, the averages are not really necessary 
but they are usually written down. In spin-glasses and glasses this observable does not 
yield information on the existence of any growing length as we shall discuss below. 

2.3 Two-time quantities 

The auto-correlation function and linear susceptibility are defined as 

NC(t,t w ) = EiLi(<?i(t,0) = Eti(Ssi(t)s Si (t w )) , 

N X (t,t w ) ee J2ti(Xi(t,t w )) = ElJi dt' R(t,t>) , 1 ) 

with Ri(t,t') = 8(si(t))/Shi(t')\ h=0 and NR = 5Z i=1 Ri the local and global instanta- 
neous linear responses, respectively. 

A number of fully general relations between the linear response and a correla- 
tion computed as an average over unperturbed system trajectories have been recently 
derived. Quite generally they read 



R(t,tw) - 



dC{t,t v 
dt w 



N 



# -1 l>i(*)Si(*.»)> 



0(t -t w ). (3) 



The explicit form of the factor Bi in the second term in the right-hand-side (r.h.s.) de- 
pends on the microscopic dynamics; it has been computed for Langcvin processes (|16[) 
and Markov processes for discrete variables (fTT)) . In all cases Bi is the deterministic 
drift in the sense that (s^ = (Bi). 

In equilibrium, the second term in the r.h.s. of Eq. ([3]) is equal to the first one, and 
R (or x) and C are related by the fluctuation-dissipation theorem 

R(t,t w ) = ±^^6(t-t w ), (4) 

X(t, t w ) = i [C(t, t) - C(t, t w )]9{t - t w ) , (5) 

and all two-time quantities depend on t — t w only. 
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2.4 Order parameter 

The very concept of a glass order parameter is not obvious. Although some kind of 
static amorphous order may develop with time in a glassy regime, searches have given 
negative results so far (not surprisingly since one does not really know what is looking 
for) (|18p . The simplest possibility, (si(t)), is void of information and the space-time 
spin-spin correlation C(r, t) (the Fourier transform of the structure factor) is, to a 
first approximation, time independent and very similar to the one in the super-cooled 
liquid ®. 

Edwards-Anderson order parameter. A dynamic order parameter, named after Ed- 
wards and Anderson (EA) who introduced it for spin-glasses, is defined as 

q ea = lim lim C(t,t w ) , (6) 

t— J-oo t w — >oo 

where C(t, t w ) is the auto-corrclation function. The order of the two long time limits is 
crucial since the weak long-term memory (2; 3) ensures that limt u ,_ i . 0O lim^oo C(t, t w ) — 
(see Fig. 11.11 (a)). The EA order parameter detects ergodicity breaking: in the liquid 
(paramagnetic) phase q ea = because the system looses memory at finite time scales 
but in the glass phase such memory remains. The idea is quite generic and it can be 
applied to glassy systems made of constituents of any kind. 

Replica overlap - a static counterpart. The theory of spin-glasses (fll?)) suggests the 
definition of a static order parameter as the overlap between two replicas (two sys- 
tems with the same quenched randomness) a and b subjected to a fictitious attractive 
coupling of strength e that forces them to be in the same thermodynamic state: 



1 N 

vE<**>«- (7) 



q = lim lim 

e-s-0 AT-j-oo N 

The spin-glass susceptibility is its linear response to an infinitesimal variation of the 
coupling: 



dq 

Xsg 7T~ 

oe 



= |E(^k 6 )-W)K 6 4)) • (8) 

£=0 iV 



The averages (. . . ) e and (...) arc taken here with the corresponding Gibbs-Boltzmann 
distributions. A crucial problem with these definitions is that below T g aging persists at 
any finite (with respect to some function of N) time scale. Nonetheless, \sg motivates 
the definition of a dynamic analogue that could detect the growth of static order during 
aging, see Sect. 12.61 

2.5 Dynamically correlated volume 

A more natural proposal is to consider the spatial variation of the speed of relaxation 
in different regions of space. This point of view is very close in spirit to the idea of 
dynamical heterogeneity in the supercooled liquid regime as characterized by local 
relaxation times. As discussed later, this point of view also emerges naturally in the 
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theoretical analysis of soft-modes associated with the time reparametrization invari- 
ance in the effective dynamical equations of motion at long time scales ©. We present 
below different ways to access a correlation related to the local dynamics. 

Spatial correlation of local dynamics. Let us regard the local two-time dependent 
auto-correlation function Ci{t,t w ) as the local order parameter. Its spatial variation 
can be easily characterized by the spatial correlation function ([T1)l [201 l2Tj) 

Gi(i,j;t,t w ) = (Cifa^Cjfatu)) - [C{t,t w )} 2 . (9) 

The subscript 4 is to remind one that it is a 4-point function since the auto-correlation 
function is already a two-point correlation function (in time). 

The key question is whether there is a characteristic dynamical length-scale, ^(t, t w ), 
beyond which G4 de-correlates, e.g. 

-lnG 4 {i,j;t,t w ) ~ /* +&lnr (10) 

?4(t, tw ) 

with r = \ fi —fj\- [This expression assumes the usual decay G4 ~ r~ b exp(— but 
more general forms of the type G4 — r~ b have also been considered.] 
By integrating Eq. (|10j) over space one obtains 

Nx4(t,t w ) = Y^Gi(i,j;t,t w ) = (§2Ci(t,t w )} 2 ) - (Y,C<(t,t w )) 2 (11) 

i,j i i 

which measures the dynamically correlated volume. Note that in spite of being usually 
called a 'susceptibility', strictly speaking X4 is not one ([221 I25| . The same proviso 
applies to Xsg in Eq. (|13p . 

The 4-point correlation function in Eq. ([9]) and the integrated one in Eq. (fTT]) are 
extensions of similar expressions defined for the supercooled (stationary) liquid state. 
In the glassy regime one needs to keep the t w dependence in ^4(t,t w ). Still, it can 
be measured in numerical simulations and in some experiments from real-time real- 
space images. One has to keep in mind, though, that this quantity detects how similar 
motion in different regions is but not necessarily whether these are ordering in the 
same state. 

Higher order susceptibilities. Although G4 has been studied in numerical simula- 
tions ([2T1 124")) its direct experimental investigation remains a challenge as, in general, 
multi-point correlators. (Lucky exceptions are colloidal suspensions in which confo- 
cal microscopy allows one to store the full particle configuration (|8]).) A natural way 
out would be to measure responses to external perturbations, actual susceptibilities, 
as suggested in ([251 12S)) and done experimentally in ([27)1. The basic idea is that G4 
should be related to a non-linear susceptibility by some sort of generalization of the 
FDT, much in the same way as the ordinary correlation function is linked to the linear 
susceptibility, see Eq. ([5j. One could then measure the latter to extract information 
on the former. In order to fulfill this program it is necessary to establish which are 
the non-linear susceptibilities associated to multi-point correlators and which is the 
generalization of the FDT, holding possibly out of equilibrium. Exact general relations 
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between multi-point correlators and non-linear response functions (|22| |23( 126)) derived 
for systems subjected to a Markovian dynamics show that beyond linear order the 
susceptibilities are related not only to multi-spin correlations (such as G4) but also to 
more complicated correlators. For instance, for the second order response of two spins 
Rtf(t,ti,te) = S 2 (s l (t)s J (t))/Sh l (t 1 )h j (t 2 )\h=o in equilibrium one has 



d_d_ 



'Si(t)Sj(t)Si(ti)Sj(t2)) 



(12) 



for ti 7^ t%. This feature poses the problem of choosing the best suited non-linear 
susceptibility to detect cooperative effects. In a series of papers both third (f2"6"]l and 
second (|23[) order susceptibilities have been considered. Analytical and numerical stud- 
ies show that the non-linear susceptibilities and G4 obey analogous scaling forms from 
which one can extract a cooperative length. Experimental studies are on the way (|28j) . 

The second order susceptibility x^HMto) = ft dti / t * dti K?' (t,ti,t2) is strictly 
related to the fluctuations of t w ) (|2"51 |2T))) . 

Distributions of coarse-grained two-time quantities. Another way of extracting a 
growing length, £, alternative to that expressed by Eq. (fT0|) . is to study the full distri- 
bution of local two-time functions (fl~5 l l30 l [3T | l32|) . Indeed, for I 2> £ the coarse-graining 
boxes naturally become independent and one should recover a Gaussian distribution. 
The crossover from non-trivial to trivial dependence can then be used to estimate £ 
that should, presumably, behave as £4. The same argument can be applied to the pdf 
of the local linear susceptibility Xi s - 

2.6 Growth of underlying static glass order? 

An intriguing problem is whether any static glass order develops during aging. The 
spin-glass susceptibility, Eq. §5$ , suggests to define a 4-point correlation function 



Nxs 9 W = J2 G ^r,t) = £< S ?(^ a (^(i)4W> • ( 13 ) 

Time t is measured after the temperature quench at which the two replicas, labeled a 
and 6, are prepared in independent random initial configurations. One is interested in 
finding a dynamic length-scale, £, S g(t), beyond which G sg decorrelates, 

-\nG sg (i,j;t)~ — T —r+clnr. (14) 

£ S g(t) is simpler than £i(t,t w ) in that it is a one-time quantity like the domain size 
L(t) in usual phase ordering processes. 

In the context of spin-glasses the value of the parameter c is used to distin- 
guish between a disguised ferromagnet picture, as proposed in the droplet model ([1]), 
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and a more complex equilibrium structure, such as the one predicted by mean-held 
models (|19p . In the former lim,,-^ G sg (r) — >■ const and c = while in the latter 
lirtv^oo G sg (r) — > and c > 0. 

It is interesting to compare Xsgif) and Xi (t,t w )- Since there are no interactions 
between the two replicas a and b, purely dynamic correlation cannot exist between 
them. Thus, the dynamical SG susceptibility Xsg{t) can only detect growth of (if any) 
static order much as domains in usual phase ordering. On the other hand X4(A t w ) can 
detect both static and dynamic order. 

2.7 Effective temperature 

Although we shall not develop the T e ff ideas here we include a short paragraph re- 
calling its definition; the concept will appear in a number of places later in the article. 

The deviation from the fluctuation-dissipation theorem (FDT) found in mean- 
field glassy models (|3l U} and in a number of finite dimensional systems (|37|) can be 
rationalized in terms of the generation of an effective temperature (|38| in the system. 
The identification of a temperature out of equilibrium makes sense in the asymptotic 
limit of small entropy production (cither long times or very weak applied drive, see 
Sec. 15. 4p in which the system evolves slowly. T e ff is defined as minus the inverse slope 
of the asymptotic parametric plot x(C), where x is the linear integrated response 
and C is the two-time correlation. The thermodynamic character of T e // has been 
checked in mean-field models and mode-coupling theories as well as in simulations 
of Lennard- Jones mixtures, models of silica and many others. A basic condition is 
that T e ff obtained from different (interacting) observables should be equal whenever 
measured in the same dynamic regime. A different scenario is found at the lower critical 
dimension and in critical dynamics, and in trap models with unbounded trap depths. 



3 Phase ordering 

Phase-ordering kinetics is an important problem for material science but also for our 
generic understanding of pattern formation in noncquilibrium systems. 

Let us consider a physical macroscopic system in contact with an external reservoir 
in equilibrium. Imagine now that one changes a parameter instantaneously in such 
a way that the system is taken from a disordered phase to an ordered one in its 
(equilibrium) phase diagram. Two paradigmatic examples are spinodal decomposition, 
i.e. the process whereby a mixture of two or more substances separate into distinct 
regions with different concentrations, and magnetic domain growth in ferromagnetic 
materials quenched below the Curie temperature ([5]). 

Closely related to the above is the process whereby a critical state is approached 
via a quench from the disorder state right to the temperature at which the phase 
transition occurs. 

In both sub-critical and critical coarsening the dynamical process starts from an 
equilibrium high temperature disordered state and progressively evolves building a 
new phase, stable at a lower temperature, either ordered or critical. 



Phase ordering 11 



3.1 Growing length, aging and scaling 

The evolution of an initial condition that is not correlated with the final equilibrium 
state (and with no bias fields) does not reach equilibrium in finite times. More ex- 
plicitly, domains of all the phases of the equilibrium state at the final temperature 
T keep on growing, until their typical size Lit) becomes of the order of the system 
size L. For any shorter time the system is out of equilibrium and, in particular, the 
non-equilibrium evolution does not come to an end whenever the thermodynamic limit 
L — > oo is taken at the outset. 

The very existence of a growing length L(t) is at the heart of the aging behavior 
observed in these systems as can be easily understood. Since L(t w ) increases with 
t w the system needs larger rearrangements to decorrelate from older configurations. 
This simple fact is the origin of the strong t^-dcpcndcncc of the autocorrelation and 
linear-response described in Sec. [T] and depicted in Fig. 11.11 

In the asymptotic time domain, when L(t) has grown much larger than any micro- 
scopic length in the system, a dynamic scaling symmetry sets in, similarly to the usual 
scaling symmetry observed in equilibrium critical phenomena. According to this hy- 
pothesis, the growth of L(t) is the only relevant process and the whole time-dependence 
enters only through L(t). Observables such as correlation and response functions take 
precise scaling forms that will be discussed in Sees. 13.21 and 13.31 Exceptional cases 
where dynamic scaling is not observed will be shortly discussed in Sec. 13.31 

3.2 Aging at a critical point 

The scaling behavior of binary systems quenched to the critical point is quite well 
understood since this issue can be addressed via scaling arguments (39) and rcnormal- 
ization group approaches (|40|) which give explicit expressions for many of the quantities 
of interest up to two loops order. Numerical simulations (|4Tj) confirm the analytic re- 
sults and probe exponents and scaling functions beyond the available perturbative 
orders. In this case the system builds correlated critical clusters with fractal dimen- 
sion D = (d+2 — i]) /2, where r\ is the usual static critical exponent, in regions growing 
algebraically as L(t) ~ i 1 / Ze «, z eq being the dynamic equilibrium critical exponent 
relating times and lengths. 

In the asymptotic time domain the correlation function ([TJ has the scaling form 



The pre-factor L{t)~ 2 ( d ~ D ^ takes into account that the growing domains have a fractal 
nature (hence their density decreases as their size grows) and the dependence on r/L(t) 
in f(x) expresses the similarity of configurations at different times once lengths are 
measured in units of L(t). 

For two-time quantities, when t w is sufficiently large one has 




(15) 




(16) 
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Here C st (r) ~ Z^r)" 2 ^- ) and i? sf (r) ~ ^(t)" 2 ^- )-^, where t = f - t„. The 
scaling functions f c and / r describe the non-equilibrium behavior and take the lim- 
iting values / c (0) = f r (0) = 1 and / c (oo) = / r (oo) = 0. The correlation and re- 
sponse function of the equilibrium state at T c obey FDT. i? s t(r) = f3C s t(r). In the 
scaling forms the equilibrium and non-equilibrium contributions enter in a multiplica- 
tive structure. Non-equilibrium effects are taken into account by taking ratios be- 
tween the sizes of the correlated domains at the observation times t w and t in the 
scaling functions. Of a certain interest is the limiting fluctuation-dissipation ratio 
T/T e fj = = limt-^oo limt^-xx, = TR(t,t w )/[dC(t,t w )/dt w }, due to its universal 
character ([39)) . 

Experiments on the non-equilibrium kinetics near a critical point are reported in 
(|42|) . where the orientation fluctuations of the director of a liquid cristal are measured 
after a sudden change of the control parameter (in this case an AC voltage) from a 
value in the ordered phase to one near the critical point where the Frecdcricksz second 
order transition occurs. In this quenching procedure the initial state is ordered. Ex- 
perimental data show a behavior of two-time quantities in substantial agreement with 
the scaling pattern described in Eqs. (|16|) . as expected since the system can be de- 
scribed in terms of a time-dependent Ginzburg-Landau equation, similarly to ordinary 
magnetic systems. Interestingly enough, even the limiting fluctuation-dissipation ratio 
Xoc turns out to be in good agreement with the value found in the two-dimensional 
Ising model. 

3.3 Aging in the low-temperature phase 

The late stage of phase-ordering in binary systems is characterized by a patchwork 
of large domains the interior of which is basically thermalized in one of the two equi- 
librium phases while their boundaries are slowly moving producing the power-law 
L(t) ~ t x / z . This picture suggests the splitting of the degrees of freedom (spins) into 
two categories, providing statistically independent contributions to observables such as 
correlation or response functions. More precisely, a quasi-equilibrium stationary con- 
tribution arises as due to bulk spins, while boundaries account for the non-equilibrium 
part. Then (j43|) . asymptotically one has 

C(r,t)~C st (r) + C ag (r,t) (17) 

The first term describes the equilibrium fluctuations in the low temperature broken 
symmetry pure states 

C st (r) = (1 - ( Si )l q ) g Q-) , (18) 

where (si) eq is the equilibrium expectation value of the local spin, and g(x) is a function 
with the limiting values g(0) = 1, linxj-^oo g(x) = 0. The second term takes into 
account the motion of the domain walls through 

C ag (r,t) = ( Si ) 2 eg f (-^) , (19) 

with /(l) = 1 and lim a: _j. 00 f(x) = 0. Both C s t and C ag obey (separately) scaling forms 
with respect to the equilibrium and the non-equilibrium lengths £, L(t). In particular, 
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Eq. (fl~9| expresses the fact that system configurations at different times are statistically 
similar provided that lengths are measured in units of L(t), namely the very essence 
of dynamical scaling. 

An analogous additive separation holds for two time quantities 

C{t,t w )~C at {T)+C ag {t,t w ), (20) 

and similarly for the response function. The equilibrium character of the first term 
implies 

C st (T) = (l-{ Si )l q )g c (^±y (21) 

where g c {x) is a function with the limiting values p c (0) = 1, lim. J ;_ s . 00 g c {x) = 0. The 
non-equilibrium term obeys 

C ag (t,t w ) = ( Si )t q f c (J^y (22) 

with / c (l) = 1 and lim a; _i. 00 f c (x) = 0. In the long t w limit the two terms in Eq. (|2"0"]) 
vary in completely different two-time scales. The first one changes when the second 
one is fixed to q ea = (si)l q , see Eq. (j6]), at times such that L(t) / L(t w ) ~ 1. In this 
regime C decays to a plateau at q ea = (s^) 2 , see Fig. 11.11 The second one varies when 
the first one has already decayed to zero. The mere existence of the second term is 
the essence of the aging phenomenon below T c with older systems (longer t w ) hav- 
ing a slower relaxation than younger ones (shorter t w ). Such a sharp separation of 
timescales is the hallmark of quenches in the ordered phase, at variance with crit- 
ical quenches where equilibrium and non-equilibrium contributions are entangled in 
the multiplicative form (|16j) . Although we have used the terminology of binary sys- 
tems, the scaling structure discussed insofar is believed to hold quite generally 
(some exceptions will be discussed below), including systems with more than two low 
temperature equilibrium phases (as described, e.g. by Potts or clock models) or with 
a continuous symmetry (i.e. vector O(N) models). Contrarily to the case of critical 
quenches, however, a systematic expansion method in sub-critical quenches is much 
more difficult (|4"4"j) and general results comparable to those at T c are not yet available. 
The scaling structure discussed above has been proven analytically only in special 
cases, such as the one dimensional Ising chain with Glauber dynamics (1391 145[) or the 
Langevin dynamics of the <i-dimcnsional O(N) model with non-conserved order pa- 
rameter in the large N limit (|46|) . It is supported by semi-analytical arguments in two 
dimensions (|47 { l48 j) . numerical simulations (|50| l5Tj) and approximate theories (|53|) . It 
has also been confirmed in experiments (|54p. 

The growing length. The growing length depends on a few characteristics of the or- 
dering process - the dimension of the order parameter, whether there are conservation 
laws, the presence of quenched disorder - and may serve to classify systems in classes 
akin to universality ones ((5]). Although the following results are hard to obtain with 
rigurous arguments, they are by now well established. Basically, one distinguishes two 
important cases: in clean models L{t) ~ t x / z ; in problems with quenched randomness, 
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where topological defects are pinned, the growing length is expected to slow down 
from a power law to a logarithmic dependence on time (due to thermal activation 
above barriers with a power-law distribution) see Sect. |4j The dynamic exponent z 
depends on the type of microscopic dynamics and dimension of the order parameter. 
For example, in curvature driven growth with scalar order parameter z = 1/2 and 
in phase separation with scalar order parameter (and no hydrodynamics) z = 1/3. A 
crossover at a static length explains how the activated scaling can be confused with a 
temperature-dependent power law in dirty cases ([55]) as explained in Sect. [4] 

The scaling functions. The full theoretical description of a coarsening process neces- 
sitates the determination of the scaling functions. This is a hard task and there is no 
powerful and systematic method to attack this problem yet. 

Still, Fisher and Huse proposed that the scaling functions in the aging regime, e.g. 
f c , should be robust ([I]). Changes in the model definition that do not modify the nature 
of the equilibrium initial nor target ordered state - such as weak quenched disorder 
not leading to frustration - should not alter the scaling functions. This is the so-called 
super-universality hypothesis. In this way the scaling functions in spin models with 
random ferromagnetic bonds or random fields should be identical to those found in 
the pure limit. Numerical tests in d > 1 systems point in the direction of validating 
this hypothesis (j48| |49| l5T|) while very recent studies of the scaling properties of the 
linear response in the d = 1 random bond ferromagnet tend to falsify it (|52p . 

The four-point correlation G4 . The behavior of G4, and the way L(t) is encoded in 
it, can be easily understood in coarsening systems. In Fig. 13.31 a configuration of the 
system with two interfaces (denoted 1 and 2, continuous lines) at time t w is sketched. 
The dashed interface denoted as 3 is the location of interface 2 at the later time t (such 
that t — t w < t w ). For t — t w <C t w one has Ci(t,t w ) = 1 - for concreteness we think 
in terms of Ising spins - everywhere except in the regions spanned by an interface in 
the interval (t w ,t) - the region between 2 and 3 in Fig. 13.31 - where Ci(t,t w ) = — 1. 
Hence G4 decays on a volume of order L d {t) — L d (t w ) which grows in time, and £4 
increases. For longer time differences (t — t w ~ t w ), however, another interface present 
at t w (1 in the sketch) may superseed at time t the position of interface 2 at the 
previous time t w . From this time on, \4 stops growing and saturates to a value of 
order L d (t w ), namely the typical volume contained between two interfaces at t w (i.e. 
between 1 and 2). The saturation can also be explained by observing that for long 
time differences, t — t w ^> t w , xa factorizes as (l/N)J2 i j(si(t)sj(t)){si(t w )sj(t w )). 
Enforcing scaling, (si(t)sj(t)) ~ f[r/L(t)], one recovers lim^^^oo Xii^i tw) °c L d (t w ). 
In this way, L{t w ) can be extracted from the long t — t w behavior of \4 but, of course, 
this is not really necessary in such coarsening problems. Importantly enough, L[t) is 
a measure of growth of static order in these cases. 

Exceptional cases: breakdown of dynamic scaling. There are also a number of coars- 
ening systems in which dynamic scaling breaks down ([5]). Among these are the XY 
model [0(2) symmetry] in one (|56|) and two (|57|) dimensions, the clock model with 
p > 4 states (58), the one-dimensional Hcisenberg model [0(3) symmetry] (|59|) and 
the large- N model with conserved dynamics (|60|) . Interestingly, in most cases this is 
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Fig. 3.3 Schematic representation of the mechanism of growth and saturation of Xt(t,t w ) 
in coarsening systems. 



due to the presence of more than one length growing macroscopically during phase 
ordering. Let us consider the XY model in d = 1 as a simple paradigm. The order 
parameter in this case is a planar vector. For long times after a quench from a high 
to zero temperature, when the excess energy is greatly reduced, the continuous nature 
of the order parameter allows only soft (Goldstone) modes, namely smooth rotations 
of the order parameter, called textures. A winding length L w can be dchned as the 
typical distance over which a 2tt rotation of the spin is detected. In addition, due to 
the symmetry of the disordered initial state and of the Hamiltonian between clockwise 
(textures) and counterclockwise (anti-textures) rotations, both topological defects are 
formed in the evolution, separated by winding inversions. Denoting with L the typical 
distance between winding inversions, it is clear that the reduction of energy implies 
the growth of both L w and L as time elapses. It was shown ([561 that these lengths 
grow with different exponents and this phenomenon prevents dynamic scaling to set 
in. 

3.4 Elastic manifolds 

The dynamics of a directed <i-dimensional elastic manifold embedded in an N dimen- 
sional transverse space is simply modelized by a Gaussian scalar field theory that 
goes under the name of Edwards- Wilkinson (EW) model. Non-linear effects can be 
accounted for by including a Kardar-Parisi-Zhang (KPZ) term. It has been known 
for some time that these systems age and that the aging dynamics is governed by a 
growing correlation length (|5} [HI [61] |62| . 

The aging regime, before saturation at a time-difference that depends on the length 
of the line, L, is characterized by a multiplicative scaling of two-time quantities. For 
instance, a two-time generalization of the roug hness, (w 2 )(t,t w ) = L d J d d x([Sh(x,t) - 
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Sh(x, t w )] 2 ), with Sh(x,t) = h(x,t) — (h(x,t)}, h the height of the manifold and x the 
position on the d-dimcnsional substrate, scales as 

(w 2 )(t,t w )^L 2 <(t w ) f w2 (j^j) ■ (23) 

Assuming lim^^co f w i (x) = x 2 ^ , in the limit r 3> t w this observable reaches a sta- 
tionary regime in which (w 2 )(t 7 t w ) ~ L 2< =(t). For even longer time-delays such that 
L(t) — > L one finds saturation at L 2< >. The roughness exponent £ is due to thermal 
fluctuations and it is simply equal to (2 — d)/2 in the EW manifold. The cross-over 
to saturation is then described by an extension of the Family- Vicsek scaling that 
takes into account the out of equilibrium relaxation and aging effects. Correlation and 
response functions are related by a modified fluctuation-dissipation relation after re- 
moving the diffusive factor and a well-defined effective temperature (|3"5|) exists and 
depends on the initial (To) and final (T) values of the temperature before and after 
the quench with T eff > T if T > T and T eff < T if T < T (ED [52]). 

4 Role of activation: the droplet theory 

The dynamics of glassy systems at low enough temperatures should be dominated 
by thermal activation. Although it is very difficult to study activated processes from 
first principles, several phcnomcnological proposals for models with quenched disorder 
exist. In particular, a droplet picture has been put forward for spin-glasses ([T] 1631 
I64p and related systems including elastic manifolds in random media (I65p and vortex 
glasses (|6"6")l . This model assumes that a static low temperature phase, associated with 
a zero temperature glassy fixed point in a rcnormalization group sense, exists in these 
systems. Droplet-like low-energy excitations of various sizes L on top of the ground 
state render the dynamics strongly heterogeneous both in space and time. The typical 
free-energy gap of a droplet with respect to the ground state and the free-energy barrier 
to nucleate a droplet are assumed to scale as L and L^, respectively, with 9 and i/j 
two non-trivial exponents. Static order is assumed to grow as in standard coarsening 
systems. Dynamical obscrvablcs such as the two-time auto-correlation function should 
then follow universal scaling laws in terms of a growing length L(t) originated in 
Arrhenius activation over barriers growing as a power of the length 

i-roe^ => m-lTHt/ro)} 1 ^ . (24) 

The strong-disorder renormalization approach in configuration space (67) yields an 
explicit construction in favor of the droplet logarithmic scaling. 

In the case of spin-glasses the very existence of static glassy states is accepted but 
their detailed nature is a much debated issue. Moreover, it is far from obvious whether 
coarsening of only two competing states occurs and whether the growth is determined 
by thermal activation over such barriers. In the present section we discuss a recipe to 
examine the droplet picture quantitatively. 

4.1 Efficient strategy for data analysis 

In practice, the asymptotic dynamic scaling features associated with the putative 
zero temperature glassy fixed point are difficult to access in numerical simulations 
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and experiments. The following strategy helps avoiding the control of pre-asymptotic 
effects that last for very long (|55|) : 

I) Measure the dynamical length L(t) and analyze it over the full time duration by 
taking care of the crossover from the initial non-activated dynamics, that could 
be diffusive, critical or else, to the asymptotic activated regime. 

II) Reparametrize the time-dependent quantity of interest, say A(t), using L(t) ob- 
tained in I) as a time-length dictionary. 

In this way pre-asymptotic corrections are dealt with separately: those due to L 
and those due to the scaling functions. Once the 'dictionary' L(t) is determined, step- 
II) is very much straightforward: no uncontrolled fits are needed since the essential 
exponents (such as the energy exponent 9) are provided by independent studies of 
static properties. We prove the efficiency of this strategy by analysing numerical results 
for Sinai diffusion (|T0[) and we recall the study of the random manifold and Edwards- 
Anderson spin-glass along these lines. 

Sinai model: a test case. The Sinai model is a random walker hopping on a one- 
dimensional lattice i = 1,2, ...N over which a quenched random potential U is de- 
fined. The statistics of the random potential are such that [(Ui — Uj) 2 ] = r^j where 
rij is the distance between sites i and j and [■ ■ • ] stands for the average over dif- 
ferent realizations of the random potential. The statistical analysis of disorder yields 
= ip = 1/2. 

The walker starts from a randomly chosen initial point at time t = 0. The mean- 
squared displacement is H(t,t w ) = [((x(t) — x(t w )) 2 )] with x(t) the position at time 
t and (• • ■ ) an average over thermal histories and initial conditions. The rigorous 
analysis of \/B(t, 0) = L(t) gives the diffusion law L(t) oc [T\n(t/ro)] 2 at temperature 
T (|70|). Also of interest is the linear-susceptibility, %(£, t w ), of the averaged position 
with respect to a small bias field acting on the walker after the waiting time t w . In 
equilibrium the FDT relates the two quantities as x(i, t w ) = B(t, t w )/2T. 

The asymptotic scaling properties have been almost fully uncovered by a real- 
space renormalization group (RSRG) approach (|7T|) which is believed to become exact 
asymptotically. For illustrative purposes, we reproduce some of these results with 
droplet scaling arguments. 

In the quasi-equilibrium regime L(t = t — t w ) < L(t w ) the RSRG predicts B(r + 
tw,t w ) = 2Tx(T+t w , t w ) ~ TL 3 / 2 (T)g b (L(r)/ L(t w )) whereof, is a scaling function. The 
exponent 3/2 can be explained as follow. In first approximation at time t w the particle 
sits in the lowest energy minimum within a length scale L(t w ) around the initial 
point. Since L(t) <C L(t w ) the particle does not have time to explore regions which 
are far from this minimum and diffusion typically vanishes in the quasi-equilibrium 
regime. However, with a small probability 8Uh/L 9 (t) ~ T/L°(t) with 9 = 1/2 there 
is a secondary energy minimum within the length scale L(t) with the energy barrier 
between them being lower than the thermal energy T. In such a rare event the particle 
can hop to the secondary minimum yielding the disorder-averaged behaviour B(r) ~ 
T j L x I 2 (t)L 2 (t) ~ TL 3 / 2 (r). 

In the aging regime L(t) > L(t w ) the particle diffuses over longer and longer lengths 
looking for lower and lower energy minima. Within a given time scale t the particle 
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moves to the right or to the left by an amount L(t). However, the direction of motion 
is almost deterministically given by the direction with the lowest energy barrier since 
the difference between energy barriers (not minima) is of order L 1//2 (i), a diverging 
quantity in the long times limit. One then has B(t,t w ) <~ L 2 (t)fi,(L(t)/L(t w )). This 
also implies that the linear susceptibility typically vanishes in the aging regime since 
a change of the potential of order SUh = hL{t) (an analogue of the Zceman energy in 
magnets) induced by an infinitesimal field h docs not affect the difference in energy 
barriers of order L x / 2 {t). However, there are rare events such that the energy barriers 
to the left and right are almost degenerate with only a small difference in height SUh ■ 
The probability to find such degenerate barriers scales as SUh /It with ip = 1/2. The 
infinitesimal h then induces a change in the direction of diffusion resulting in a disorder- 
averaged displacement of order (SUh/L^)L and thus x{% tw) ~ L 3 / 2 (t)f x (L(t)/L(t w )). 




t log (t/-uO(T)) 



Fig. 4.4 Diffusion in the Sinai model: the mean-squared displacement from MC simulations, 
a) Raw data at T — 1.8, 1.6, . . . , 0.4, 0.2 from top to bottom. The unit of time t is one MCs. 
The average is taken over 10 5 realizations of the random potential, b) Scaling plot of the 
dynamic length L 2 (t) = B(t, 0) showing the crossover behaviour from normal diffusion at 
short time t - see the inset - to the expected L 2 ~ ln 4 (f/ro) asymptotic law. 

We performed Monte Carlo (MC) simulations to examine the anticipated asymp- 
totic scaling IrJ^)) . In Fig. 14.41 a) we show the mean-squared displacement B(t, 0) 
of the Sinai walker at different temperatures, decreasing from top to bottom. It de- 
fines unambiguously a dynamical length L(t) = \/B(t, 0) that we shall use below to 
reparametrize time dependent quantities. The data demonstrate that extremely long 
times are needed to reach the L(t) ~ In t/ro asymptotic result known to be exact 
analytically. At finite T the short-time diffusion is normal B(i,0) = Dt with a diffu- 
sion constant D, see the inset in Fig. 14.51 The crossover to the activated regime takes 
place at a 'thermal length scale' L (T) such that SU(L (T)) = L^ 6 (T) ~ T and the 
corresponding time scale tq(T) is fixed using Lq(T) = Dtq(T). As shown in Fig. 14.41 
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Fig. 4.5 Diffusion in the Sinai model: two-time mean-squared displacement and linear sus- 
ceptibility from MC simulations at T — 0.6. Note the absence of a developing plateau, a) Raw 
data of the mean-squared displacement B(t,t w ) (open symbols) and the linear-susceptibility 
x(t,t m ) (filled symbols), b) Scaling plots in terms of the dynamic length. Inset: zoom over 
the quasi-equilibrium regime. 

b) the slow crossover is well described by proposing that the scaled length L{t) / Lq(T) 
is a universal function of the scaled time £/to(T). The asymptotic Sinai's diffusion law 
can be parametrized precisely with L(t) ~ Lo(T)[\n(t /tq(T))] 2 . 

Next we proceed to step II) . In Fig. 14.51 a) we show the aging effects observed in 
the mean-squared displacement and the corresponding linear-susceptibility. In Fig. 14.51 
b), the numerical data are reparamerized by the dynamical length scale obtained in 
the step I) and follow the expected scaling laws. This means that the correction to the 
asymptotic behaviour is not large in terms of length scales. The absence of a plateau in 
B and x m their raw and scaled forms demonstrates the fact that there is no additive 
separation of time-scales - as in Eq. (|20|) - in this problem. 

Elastic manifolds in random media. The dynamics of a directed d-dimensional elas- 
tic manifold embedded in an N dimensional transverse space under the effect of a 
quenched random potential play an important role in a variety of physical systems 
ranging from coarsening in dirty systems to fracture. An application is the one in 
which the directed lines are vortices in super-conductors aligned in the direction of 
the magnetic field and simultaneously pinned by impurities. 

These systems are intimately related to the problem of Sinai's diffusion. The sim- 
plest case is a 1+1 dimensional directed polymer in random media (DPRM). Each 
configuration is described by the transverse displacement x(z, t) at position z along 
the directed polymer at time t. Two natural quantities used to characterize aging are 
the mean-squared displacement B(z, i, t w ) = [((x(z, t) — x(0, t w )) 2 )} and the linear sus- 
ceptibility x( z ,t,t w ). The FDT x( z ,t,t w ) = 3(z,t,t w )/2T is satisfied in equilibrium. 

Scaling arguments of droplet type naturally apply to the present case (|68|). To this 
end one just needs to keep in mind the scaling relation between the longitudinal length 
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L and transverse length L± oc L 1 * with £ the roughness exponent. Depending on the 
explored length scale, the latter takes a thermal, £t = 1/2 or a disorder, = 2/3, 
value with the crossover given at a static temperature dependent length Lq(T). The 
disorder-dominated roughness exponent is related to the energy exponent 9 by an 
exact scaling relation 9 = 2(d — 1 implying 9 = 1/3. It is also conjectured that ip = 9. 
The time evolution of a local segment x(z, t) may be viewed as diffusion of a particle 
which feels an effective (renormalized) potential created by the rest of the system, the 
variation of which scales with the transverse length as AU ~ L e ~ L°l^ D ~ h\ 1 ' t * D 
in the disorder dominated regime (|72|) . 

Using the analogy with the Sinai model one derives the scaling behaviour of dif- 
ferent observables. At a time t after the quench, the system is equilibrated up to 
longitudinal length L(t) over which the energy is higher than the equlibrium one by 
an amount of order L . Thus the energy density per unit longitudinal length, e(t) = 
U(t)/L, is expected to decay as e(t) = e(oo) + ct/L^t) 1 ^ 6 . In the quasi-equilibrium 
regime the consideration of the degeneracy of minima suggests B(z, r + t w , t w ) = 
2Tx(z,r + t w ,t w ) = TL{T)g b {z/L{t w ),L{ T )/L{t w )) [the prefactor TL(t) is due to 
T/L 9 (t)L 2 <(t) = TL(t) since 6 = 2( - 1.] In the aging regime we find B(z,t,t w ) = 
L(t) 2 <f B {z/L(t w ),L(t)/L{t w )) while x (z, t, t w ) = L(t)f x (z/L(t w ),L(t)/L{t w )) due to 
the degeneracy of barriers. 

In order to test the above, the simplest protocol is to choose a flat initial condition, 
i.e. x{z, 0) = for all z. During isothermal aging the roughness of the system develops 
progressively from short to long wave lengths. The dynamical length L(t) can be 
extracted from the growth of the static roughness B(0,t,t) = [((x(z,t) — x(z,t) 2 )} = 
B eq (z = L(t)). Here B eq (z) = lim^oo B(z, t, t) is the equilibrium roughness, which 
can be computed numerically with a transfer matrix method in the 1+1 case (|68j) . 
Another way to determine L(t) is given in (|55|) . 

By performing Monte Carlo simulations of the 1+1 DPRM the analysis of step I) 
and II) can be done precisely ([68| l55|) helped by the knowledge of the exact values of 
the exponents (+ = 1/2, £r> = 2/3 and 9 = 1/3. Quite interestingly the variation of 

1 /2 

the energy scales as AU ~ Lj_ just as in the Sinai model discussed before. 

Concomitantly with the change in roughness exponent from thermal to disor- 
der dominated values, L(t) exhibits a gradual crossover from pure diffusion with 
L(t) ~ i 1 / 2 (and z = 2) to an activated regime consistent with an algebraic growth of 
barriers (|65[) . This is similar to what is found in the Sinai model (see Fig. l4.4j) but with 
ip = 1/3. The crossover occurs at a static temperature-dependent correlation length 
Lq(T). In the analysis of numerical data the two regimes tend to be confused into 
a single one with an effective temperature-dependent power-law, L(t) ~ t 1 / z ( T ) ([g^|) . 
the T-depcndence of which is inherited from the one in Lq(T) (55), but this is just an 
approximation. Further support to the asymptotic logarithmic scaling is given by the 
rcnormalization group study in (|67|) . 

As regards the scaling properties of the two-time quantities B(z, t, t w ) and xi z > ^ tw) 
in the quasi-equilibrium and aging regimes one also faces the difficulty of going be- 
yond the thermal regime and reaching, for sufficiently long time-scales, the disorder 
dominated one. In (|55p it was shown that in the early effective power-law regime, 
approximated by L(t) ~ t 1 / 2 ^', and using the thermal roughness exponent, two-time 
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linear response and correlation functions conform to the scalings discussed above and 
a finite effective temperature exists [basically because 2(t = 1 and the diffusive pref- 
actors in B and x are both equal to L(t)]. Simulations entering the trully activated 
regime suggest that the dynamic scaling above with the disorder roughness exponent 
(d sets in (and the effective temperature progressively vanishes) (|68|) . 

Spin glasses. The isothermal aging of the Edwards- Anderson (EA) model was ana- 
lyzed along the same steps. The model is supposed to have a finite temperature phase 
transition and thus an additive separation of time-scales of the form (|20j) - the precise 
nature of the aging term depending on a two-state droplet picture ([T] 1631 I64p or a 
more complicated dynamics of the Shcrrington-Kirkpatrick (SK) type ((4]). 

The growing length is extracted from the real- replica overlap, see cq. (Q~4|) (|25|) . 
It is common to use the fit £ S g(t) ~ t 1 / 2 ^' with z(T) an effective temperature de- 
pendent exponent ([Ml 13"3l IM1 133]) . The best currently available data are from the 
JANUS collaboration suggesting z(T) ~ z(T c )T c /T flM}. However, £ ag (t) should ex- 
hibit a gradual crossover from critical dynamics at short-times to activated dynamics 
asymptotically (H2 [73]). The crossover is expected at Lq(T) ~ £ and tq(T) ~ £ Ze<! 
with the equilibrium correlation length £ = \T — T c \~ v and z eq and f the usual criti- 
cal exponents at the critical temperature T c . The values of T c , z eq and 9 are fixed 
with equilibrium studies. In the rest of the discussion we call L(t) the growing length 
-* i(t)]. 

Similar arguments to the ones exposed for Sinai diffusion and the directed manifold 
imply that the energy density per spin should decay as e(t) = e(oo) + ct/L(t) d ~ e and 
provide scaling forms for the self-correlation and linear susceptibility of very similar 
type to the ones in previous subsections (but with different values of the exponents). 

Numerical tests of the droplet picture predictions were done by a number of 
groups. The energy density decay in the d = 3 and d = 4 EA models were checked 
in ([35jl and ([36| . respectively. The numerical spin auto-correlation function C(t,t w ) 
and the linear susceptibility x(Mt«) were analyzed along the proposed scaling forms. 
In the quasi-equilibrium regime where the FDT C(r + t w ,t w ) = T\{t + t w ,t w ) 
holds, the expected scaling C(r + t w ,t w ) = q ea + ct T / L(r) e g c (L(r) / L(t w )) with 
g c (A) a 1 - const \ d ~ + . . . was verified in d = 3 (35) and d = 4 (36) (the 
fact that L(t) is far from a logarithm was attributed to pre-asymptotics and the 
analysis was perfomed using the method sketched above.). Evidence for the valid- 
ity of droplet aging scaling for the correlation function C ag (t,t w ) = f c (L(t) / L(t w )) 
and the field cooled susceptibility in which the sample is cooled in a field that it 
switched off at t w , X fc(t, t w ) = X (t, 0) - X (t, t w ) = L(t)- f x (L(t)/ L(t w )) with x (t, 0) = 
Xd — const T/L(t) s §Q in d = 4 was given in ([36]) (note that ip ~ is assumed 
here). These results were used as support to the standard droplet theory. However, 
the asymptotic value of the dynamical susceptibility \d is significantly higher than 
Xea = — qea) ([36" [ . This is clearly at odds with the standard droplet theory ([TJ 
and suggests the existence of excessive contributions from some unknown soft-modes 
other than the droplets in equilibrium. Moreover, it was found that the FDT holds up 
to significantly longer time-differences than time-translational invariance (|36|) . These 
observations suggest that, even within a droplet perspective, modifications of the con- 
ventional approach are needed. 
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A different school of thought (f34) |24|) pushes for a picture a la mean-field (namely, 
the dynamics of the SK model (j4])) that should be accompanied by a complex relax- 
ation of correlations and susceptibilities in multiple time-scales, presumably linked to 
multiple length-scales. So far, this ultramctric organization of time-scales has not been 
found numerically (|2T1 l24l [87]) . 

4.2 Other random systems 

Let us finalize this Section by briefly recalling studies of growing lengths in two other 
systems with quenched randomness: the random ferromagnet and the sine-Gordon 
model with random phases. 

On the basis of a droplet theory with barriers increasing as a power-law of dis- 
tance, the random ferromagnet should coarsen with a logarithmically increasing typical 
length (155)) . Early Monte Carlo simulations (89) showed agreement with this predic- 
tion but more recent numerical results from Ricgcr's group were interpreted, instead, 
as evidence for a power- law with variable exponent (90). This result would imply a 
logarithmic dependence of barriers on the domain size. It was argued in (|55[) that the 
effective power law maybe the effect of a very long crossover from curvature driven to 
activation with conventional power- law growth of barriers with size. Recent simulations 
support this claim (|9ip - although they were performed in a dilute ferromagnet. 

The Cardy-Ostlund or random sine-Gordon model describes the relaxation dynam- 
ics of 2d periodic elastic manifolds under the effect of a quenched random potential. 
The peculiarity of the model is that disorder induced barriers grow only logarithmically 
with size, ~ In L. A 'super-rough' or 'marginal' glassy phase exists below a critical tem- 
perature T g . The aging dynamics were studied with Coulomb gas and renormalization 
methods close to T g , and with the functional renormalization group (FRG) (|92|) and 
numerical simulations in the full low-T phase (|93|) . A dynamic length L(t) ~ t 1 / z ( T ) 
with z(T) ~ 2 + c/T develops in time, consistently with an Arrhenius argument of 
the kind described in Sect. 2J Intriguingly, = z was also found analytically and 
numerically. 

5 Growing length-scales in aging glasses 

In this Section we summarize predictions from mean-field theory as well as the outcome 
of measurements of two-time lengths in several glassy systems. 

5.1 Mean- field models 

It has been argued that the Langevin dynamics of mean-field disordered models 
are equivalent to mode-coupling theories of glasses. More precisely, the Sherrington- 
Kirkpatrick (SK) spin-glass has a dynamic transition similar to the one found in so- 
called type A mode-coupling while p-spin models with p > 3 have a random-first-order 
transition (RFOT) similar to the one found in type B mode-coupling theories (82). The 
asymptotic out of equilibrium relaxation of these models, when the thermodynamic 
limit has been taken at the outset, approaches a region of phase space that is not the 
one where equilibrium configurations lie. In models of the RFOT class this region was 
named the threshold since its energy density is at 0(1) distance from the equilibrium 
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Fig. 5.6 The two-time dependent growing length in a colloidal suspension. The data is taken 
from ((8}. </> is the packing fraction and r m ax is a cut-off used to compute the integrals (see ((8)) 
for more details). 

one ([3]). In models of the SK or type A class, although the configurations visited dy- 
namically are not the ones of equilibrium their energy density is the same ((4]). In both 
cases the region reached dynamically is fiat, in the sense that a dynamic free-energy 
landscape can be defined and its geometric properties studied, and one finds that the 
eigenvalues of its Hessian mostly vanish (J3J. The fiat directions provide, on the one 
hand, channels of aging relaxation and, on the other, their associated zero-modes give 
rise to diverging susceptibilities. In the case of mean-field models the latter cannot be 
directly associated to diverging length-scales since these models do not contain any 
notion of distance. Having said this, if one proposes that the same mechanism, namely 
the relaxation to a fiat region of phase space, is at the root of aging phenomena in finite 
dimensional glassy systems, the diverging susceptibility should be linked to a diverging 
length scale. In the super-cooled liquid regime this argument was used by Franz and 
Parisi to argue for the existence of a diverging length scale within the RFOT scenario 
from the analysis of the p-spin model (|94|) . see the chapter by Franz and Semcrjian. 
These ideas were extended to the actual MCT in (|95|) . In the trully glassy phase a 
very similar mechanism is at work and gives support to the symmetry argument © 
that we discuss in Sect. [6] and provides a scenario for dynamic fluctuations in aging 
glassy systems. 

5.2 Measurements 

In Sect. I2.ll we already mentioned simulations (fT3|) and experiments ([14]) in which the 
tagged motion of particles was followed with the intention of characterizing clusters of 
more or less mobile ones in aging glasses. 

A two-time dependent dynamic growing length, as defined in Eq. (jlOp . was mea- 
sured with numerical simulations of aging soft spheres and Lennard- Jones mixtures (1201 
l83|) . spin- glasses ([TBI [21] l24f and with confocal microscopy data of colloidal suspen- 
sions (0 l96|) . An example is shown in Fig. 15.61 where we display the outcome of the 
data analysis performed in (j8]). 



24 



In all cases £4 is described by 



£4,(t, t w ) ~ 



£,st(t — t w ) 
C «(C) 



C>q, 
C<q t 



ea. 



ea 



with k(C) a monotonically decreasing function of C and the exponent a taking a very 
small value. It is important to stress that a clear identification of a finite Edwards- 
Anderson parameter in some of these systems is hard and the separation of time scales 
might not be as clear-cut as expressed in the formula above. 

5.3 Growth processes in kinetically constrained models 

Kinetically constrained models were originally introduced in the 80s by Fredrickson 
and Andersen (97) as toy models for understanding glassy dynamics. Generally these 
are stochastic models with Markovian dynamics obeying detailed balance with respect 
to a (usually) trivial energy function. Some constraints prevent particular local tran- 
sitions between configurations (98). Despite their trivial equilibrium measure, they 
capture many features of real glass-forming systems summarized in the chapter by 
Garrahan, Sollich and Toninclli. 

The so-called spin-facilitated Ising models, are described by the Hamiltonian H = 
— rij, where Hi = 0, 1 are spin variables on a d-dimensional lattice which can be 
regarded as a coarse grained density of moving particles. The basic idea is that the 
rearrangement in a certain region around i will be facilitated if it is surrounded by 
a relatively low-density neighbourhood. In some of these systems, the existence of a 
growing length in the non-equilibium evolution after a quench can be easily recognized. 
Some time after a sudden change in temperature towards T ~ all sites are in a dense 
state with = 1 and large frozen regions, the size of which increases as the number 
of sites with m = decreases, are present ([99]) . 

A particularly interesting kinetically constrained system is the 2d spiral model. It 
has the peculiar feature of having a glass-jamming transition at the critical density 
(rii) = n c of directed percolation (|100j) . In equilibrium above n c a finite fraction of 
sites are frozen by the kinetic constraint. This means that the configuration space is 
partitioned into mutually inacccsiblc regions; namely, the dynamics is reducible. Then 
all configurations without blocked regions are dynamically connected to the empty 
state with ni = Vi, usually denoted as the high-temperature partition, and they are 
disconnected from any partition containing frozen regions. Therefore, quenching the 
system from (say) the empty state to the jammed phase, the dynamics do not freeze 
and remain confined in the high temperature partition by means of a coarsening process 
with a growing correlation length which resembles in some sense the one observed in 
binary systems without kinetic constraints (|10ip . 

5.4 Interplay between aging and drive 

The classical and quantum dynamics of large systems under an external, non-conservative, 
drive is a subject of active research. An example in the first family is the analysis of 
the rheological properties of soft glassy materials, that are a result of a competition 
between the response to the shearing forces and their intrinsic slow dynamics. An ex- 
ample in the second class is the study of quantum magnets under a current generated 
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to the coupling to leads at different chemical potential. We shall recap here and in 
Sect. 15.61 a few issues in this field that are related to the growth of a length and how 
this is affected by the drive. 

Flow in coarsening systems. The effect of an external drive on coarsening systems is 
not completely understood and in many cases a general consensus is lacking. This is 
at odds with the wide technological interest of these and similar systems in various 
application areas (p.02). In particular, the case of binary systems in a plane shear flow 
(Couette flow) has been thoroughly studied recently. In this case domains grow elon- 
gated and stretched along the flow direction. This causes ruptures of the network (| 103[) 
that may render the segregation incomplete. The characteristic lengths in the parallel 
or perpendicular directions to the flow could keep on growing indefinitely as without 
shear or, due to domain break-up, the system may eventually enter a stationary state 
(similarly to what happens in the mean field models discussed above) characterized by 
domains with a finite thickness. In general, there is no agreement on this point neither 
from experimental (|104p nor theoretical points of view. On the theoretical side, large- 
TV calculations f|105[) show the existence of an asymptotic regime with ever growing 
lengths in all directions. Stabilization of the domain size by rupture is excluded also by 
different approaches based on approximate theories for scalar systems complemented 
with RG analysis ( |106p . although with the unexpected feature of domains growing 
in the flow direction while shrinking in the perpendicular one in d = 2. Numerical 
simulations fjlOTf) cannot still establish a clear-cut evidence due to finite size and dis- 
cretization effects. It was also argued (|108|) that hydrodynamic effects are crucial in 
stabilizing a stationary state, while in the diffusive regime domains might keep on 
growing. 

Another difference with driven mean-field models is provided by the fluctuation- 
dissipation relation x(C)i the exact computation of which in the large- N model f|109[) 
excludes a simple scenario with a single effective temperature. Much more details are 
given in the contribution by Barrat and Lemaitre. 

Flow in glassy systems. The Langcvin dynamics of mean-field disordered models 
under non-potential forces has been studied in a series of papers (|110( I53|) . In all 
cases (apart from the spherical model with two-body interactions at zero temperature) 
aging is suppressed and a stationary state is reached. Two-time observables decay in 
two steps and the relaxation time decreases with the strength of the drive as in shear 
thinning systems. In the limit of weak drive the effective temperature takes two values, 
the one of the bath in the first relaxation step and a higher value in the second slower 
stage. These results have been verified in molecular dynamic simulations of Lennard- 
Jones mixtures 

The fact that a weak non-potential perturbation can change the relaxation of a 
classical system so dramatically, by rendering the dynamics stationary, can be under- 
stood in different ways. First, the type of non-potential force used in (|110|) can be 
loosely associated to an infinite temperature noise. Second, and more importantly, the 
mechanism for fluctuations based on time-reparamctrization invariancc gives a reason 
for this fact, as discussed in Sect. [5] 

Experimentally, the interplay between aging and flow in an aqueous suspension of 
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Laponite has been studied by diffusive wave spectroscopy ()112j) and light scattering 
echo experiments ()113[) (multiple scattering) and dynamic light scattering in the single- 
scattering regime (| 1 14|) . Using the latter technique Di Leonardo et al measured the 
density auto-correlation function and found that as long as the characteristic relaxation 
time is smaller than the inverse shear rate aging in unaffected by the perturbation but, 
when the relaxation time reaches the inverse shear rate, aging is strongly reduced. 

A recent study of microscopic structural relaxation of colloidal suspensions under 
shear together with relevant references appeared in (|115|) . Numerical studies of aging ad 
plastic deformations have been performed by several groups, see the reviews in (|1161 
I117p . Nevertherless, as far as we know, an analysis of growing dynamic lengths, as 
defined from G4 or x<t has not been performed in driven systems yet. 

Aging at the depinning transition of clastic manifolds. A series of analytic and numer- 
ical papers have recently showed that a non-steady critical regime, limited only by the 
steady correlation length or the system size, exists at the zero-temperature depinning 
transition of an elastic manifold in a disordered medium (fTT8| . The FRG analysis (|92| 
as well as molecular dynamics studies prove that there is a driven transient aging 
regime that displays universal features in which a growing length L(t) ~ t 1 ^ partic- 
ipates in the critical dynamic scaling description of the dynamics. At the depinning 
transition the effective height of the barriers vanishes justifying the power-law scaling 
the effect of disorder being a reduction of 2 to a value smaller than 2 (at the existent 
one- loop FRG calculation). 

5.5 Granular matter 

The gently driven dynamics of dense granular matter - athermal systems - shows many 
points in common with the ones of glass forming liquids (|119p . Aging in a water- 
saturated granular pile submitted to discrete taps has been reported in (|120p . The 
data were obtained using multispccklc diffusive wave spectroscopy to measure particle 
displacements. Evidence for dynamic heterogeneities and a growing correlation length 
in a horizontally vibrated amorphous assembly of hard disks close to the jamming 
transition was given in (|121|) . Presumably, a two-time growing length would also exist 
in such athermal systems if studied in the aging regime. On the theoretical side, the 
aging properties of a periodically perturbed mean-field model of the RFOT type and, 
more precisely, how these depend on the amplitude and frequency of the drive were 
studied in (TT22|) . 

5.6 Quantum fluctuations 

The study of the dynamics of quantum systems has been recently boosted by the devel- 
opment of cold atom experiments and advanced measurement techniques. Quantum 
magnets are expected to undergo coarsening dynamics in the ferromagnetic phase. 
The study of a mean-field model (freely relaxing and driven by an external current) 
has revealed an extension of the super-universality hypothesis to the quantum realm 
once small scales are well separated from large ones, the former feeling the quantum 
fluctuations and the latter describing the motion of large objects - the walls - being 
basically identical to the classical ones. More precisely, in O(N) field theory in the 
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large N limit or models of rotors in interaction the length over which static order is 
established grows as Lit) ~ t x / z with z = 2 as in the classical limit. This type of 
coarsening survives the application of a voltage difference by two leads at different 
chemical potential (|123|) . This prediction is open to experimental tests. 

6 A mechanism for dynamic fluctuations 

In a series of papers Chamon et al proposed that the symmetry that captures the 
universal aging dynamics of glassy systems is the invariance of an effective dynamical 
action under uniform rcparamctrizations of the time scales ([17)1 |2"T1 1501 I3"T1 I5"2")) . see 
also (124). This approach is reviewed in ((9]). Such type of invariance was first encoun- 
tered in the mean-field equilibrium dynamics of spin-glasses (|125[ I126j) and it was later 
found in the better formulated out of equilibrium dynamic of the same models ([3]|4|). 
The invariance means that in the asymptotic regime of very long times the solution 
can be found up to a time-reparametrization transformation. Within such a scenario, 
extended to systems in finite dimensions, the remaining symmetry is responsible of 
the dynamic fluctuations. Support to this claim was given by the proof of global time 
reparametrization invariance in the long times action of short-range spin-glasses (|30( 
I15j) . In the proof one assumes that the dynamics is causal and that there is an additive 
separation between a time regime with a fast relaxation and another in which it is 
slow [just as in Fig. 11.11 and in Eq. (|20p] . The invariance can then be used to describe 
dynamic fluctuations in spin-glasses and, as conjectured, in other glassy systems as 
well. The fact that the dynamic action be symmetric under uniform, i.e. spatially in- 
dependent, rcparamctrizations of time variables [t — > h(t)] suggests that the dynamic 
fluctuations that cost little action are those associated to space dependent, long wave- 
length, reparametrizations of the form t —> h(r,t). These should be the Goldstone 
modes associated with breaking time-reparametrization invariance symmetry. 

In the slow and aging regime the two-time correlation and linear response depend 
on both times and time-translation invariance is lost. As we argued in Sect. [1] the 
two-time dependent correlation acts as an order parameter. Dynamic fluctuations are 
such that ages can fluctuate from point to point in the sample with younger and older 
pieces (lower and higher values of the correlation) coexisting at the same values of 
the two laboratory times: this has been named heterogeneous aging. By looking now 
at spatially heterogeneous rcparamctrizations, we can predict the behavior of local 
correlations and linear susceptibilities and the relations between them. Different sites 
can be retarded or advanced with respect to the global behaviour but they should all 
have the same overall type of decay. Similarly, the relation between local susceptibilities 
and their associated correlations should be identical all over the sample (fl5)) leading to 
a uniform effective temperature (38). The reason why the scaling functions should not 
fluctuate much is that these are massive modes. Within the time-reparametrization 
invariance scenario the growing correlation length is due to the approach to the long- 
time regime in which the symmetry is fully dcvelopped, and the long wavelength modes 
eventually become masslcss. 

This hypothesis was put to the test in glassy systems with Monte Carlo simulations 
of the Edwards- Anderson model (Tl5l l2"Tj) . molecular dynamics of Lennard- Jones mix- 
tures (|83| , numerical studies of kinetically constrained models (j3Tj) , and in the study 
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of solvable ferromagnetic models with the analysis of the O(N) model in the large 
N limit (|32j) and the spherical ferromagnet ([84]). A scries of stringent tests some of 
them performed in these papers were summarized in The outcome of these stud- 
ies is that while trully glassy systems conform to the consequences of the hypothesis, 
simple coarsening as developed in the 0(N) and spherical ferromagnet does not with 
time-reparametrization invariance being reduced to time rescaling at the heart of the 
difference. This result is very suggestive since it implies that the invariance properties 
and the fluctuations associated to them are intimately related to the behavior of the 
global effective temperature (finite against infinite) in the aging regime. 

In short, this idea proposes a reason for the development of large dynamic fluctu- 
ations in glassy dynamics. 

7 Closing remarks 

In aging systems, the relaxation of typical features in a sample of a certain age t w 
takes a time that increases with t w . Although the reason for this fact is not necessarily 
the growth of those features, the widespread observation of growing lengths in physical 
systems, as reviewed in this chapter, promotes this one as the basic mechanism at the 
heart of most aging phenomena. (See (|127j) and the chapter by Franz and Semerjian 
for rigurous relations between growing lengths and growing times.) Related to this 
fact, the occurrence of dynamic scaling provides a parametrization of the kinetics in 
terms of typical lengths, giving access to universal properties. 

Besides these unifying concepts, the very nature of the growing length, the growth 
mechanisms, the relation with equilibrium properties etc. are issues which, in most 
cases, deserve clarification. In particular, the interplay between the growth of static 
order - as, in coarsening systems - and purely dynamical correlations - as in most 
kinctically constrained models - could help in understanding the behavior of glasses. In 
the case of spin glasses, moreover, the comprehension of the kind of order established 
within the correlated volume is fundamental to distinguish between competing scenarii, 
namely droplet-like and a la mean field pictures. 
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